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ABSTRACT 


This  thesis  describes  the  conversion  of  a  computer 
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magnitude  of  the  refractive  index  structure  parameter, 
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I.  IHTBODOCTIOH 


Many  areas  of  modern  technology  are  concerned  with  the 
propagation  and  distortion  of  electro-magnetic  radiation  in 
the  atmosphere.  In  particular,  atmospheric  optical  turbu¬ 
lence  seriously  degrades  visible  light  astronomy,  satellite 
photography  and  laser  propagation. 

Studying  and  attempting  to  predict  atmospheric  optical 
index  of  refraction  fluctuations  has  been  going  on  for  a 
number  of  years  [Ref.  1].  With  the  ascendancy  of  modern 
computer  methods  during  the  last  twenty  years,  the 
predictions  have  grown  significantly  in  accuracy  and . commen- 
surately  in  complexity.  R.E.  Hufnagel  published  a  model 
[Ref.  2],  in  1974  based  on  empirical  data  collected  at 
astronomical  observatories.  However,  it  did  not  consider 
atmospheric  dynamics  other  than  an  average  wind  speed.  In 
1577,  T.E.  VanZandt  published  a  new  model  [Ref.  3].  This 
model  is  the  most  comprehensive  and  accurate  one  developed 
to  date  and  attempts  to  include  actual  temperature  and  wind 
shear  gradients  that  force  the  production  of  optical 
turbulence. 

A.  BACKGROUND 

The  computer  code  based  on  VanZandt’ s  model  was  written 
by  J .  Warnock  at  the  National  Oceanographic  and  Atmospheric 
Administration’s  Atmospheric  Dynamics  Aeronomy  Laboratory  in 
Boulder,  Co.  The  latest  version  of  the  program,  that  is 
constantly  being  updated,  has  a  modular  format  and  consists 
of  twenty  subroutines  and  functions.  The  code  was  written 
to  be  executed  on  NOAA’s  CDC  Cyber  750  computer  using 
Fortran  IV.  The  primary  concern  of  this  thesis  was  the 
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translation  of  the  subroutines  as  received  to  code  that 
would  be  compatible  with  the  Fortran  IV  used  at  the  Naval 
Postgraduate  School’s  (NPS)  IBM  3033  computer  system  and  to 
exercise  this  program  on  specific  atmospheric  wind  and 
temperature  profiles.  Additionally: 

1.  a  main  program  was  written  to  consolidate  the  subrou¬ 
tines,  read  the  input  data  file  and  print  the  tabular 
output, 

2.  a  data  processing  subroutine  was  written  to  transform 
the  raw  input  data  into  the  variables  required  in  the 
following  subroutines,  and 

3.  a  Disspla  routine  was  written  to  produce  a  graphical 
presentation  of  the  computer  model  output  and  to  permit 
comparison  with  whatever  reference  values  were  input. 


B.  GOALS 

With  the  recent  interest  that  the  Department  of  Defense 
has  shown  in  visible  wavelength  laser  systems,  an  accurate 
computer  model  for  atmospheric  optical  turbulence  is  essen¬ 
tial  for  meaningful  systems  analysis.  Ideally,  this  model 
should  be  relatively  fast  and  easy  to  use  and  should  be  able 
to  operate  with  readily  available  data. 

The  program  that  has  teen  developed  is  capable  of  satis¬ 
fying  these  requirements.  The  raw  input  data  comes  from 
conventional  rawinsonde  printouts.  These  are  available  in 
the  computer  record  archives  of  the  National  Weather  Service 
for  hundreds  of  stations,  ashore  and  afloat,  and  going  back 
quite  a  number  of  years.  The  program  produces  tabular 
and/or  graphical  output  from  the  rawinsonde  data  entered  in 
the  data  input  file  (hereafter  DATAIN  file).  For  a  given 
location  and  date,  the  program  estimates  the  turbulence 
structure  parameter,  Cn2,  between  each  input  data  altitude. 


Because  it  is  rather  easy  to  use,  it  is  hoped  that  the 
program  will  prove  to  be  beneficial  for  physical  meteorology 
and  laser  technology  classes.  It  is  also  very  instructive 
to  modify  the  input  data  to  determine  the  effects  of  several 
atmospheric  variables  upon  the  turbulence  structure  param¬ 
eter.  For  example,  it  is  possible  to  determine  the  effect 
of  variations  in  the  temperature  or  wind  shear  profiles  upon 
the  final  output. 

The  program  should  be  useful  as  a  reference  to 
researchers  who  are  trying  to  measure  the  Cn2  profile  exper¬ 
imentally.  Assuming  the  availability  of  a  nearby  rawinsonde 
launch  site  (within  50  to  100  miles),  the  computer  model 
will  produce  Cn2  estimates  for  comparison  with  experimental 
data.  Ideally  the  experimental  measurement  device  or  the 
model  will  be  fine-tuned  as  a  result  of  detailed  comparative 
analyses. 

Finally,  the  program  also  has  the  potential  to  be  used 
as  a  subroutine  to  estimate  the  Cn2  profile  for  other 
programs.  For  example,  the  G.U.T.S.  program  presently  used 
in  AF-4706  makes  use  of  a  rather  crude  approximation  for 
Cn2.  Using  the  much  more  realistic  values  obtained  from 
this  computer  model  would  serve  to  improve  the  accuracy  of 
the  G.U.T.S.  output. 

The  primary  goal  of  my  work  was  to  translate  the  CDC 
code  into  a  smoothly  running  and  error-free  program  for  the 
IBM  system  at  NPS.  Also,  the  format  of  the  program  as 
received,  was  somewhat  rough  and  required  a  significant 
amount  of  time  to  understand  and  be  able  to  use.  As  a 
secondary  goal  then,  it  was  desired  to  make  the  program  as 
’user  friendly’  as  possible  so  that  it  could  be  used  for 
class  projects  and  by  researchers  with  a  minimum  of  instruc¬ 
tion  and  difficulty.  Finally,  since  the  Cn2  parameter  can 
vary  erratically  with  altitude,  it  was  thought  that  simple 
tabulated  output  would  be  of  limited  utility,  especially  as 


an  instructional  aid.  Therefore  a  graphical  presentation 
capability  was  desired.  Additionally,  it  is  quite  likely 
that  a  researcher  using  the  program  would  find  it  useful  to 
be  able  to  compare  several  profiles,  perhaps  obtained  by 
different  methods.  Thus  the  capability  to  read  several  data 
files  and  plot  the  comparative  data  was  also  incorporated  in 
the  program. 


II.  APPROACH  TO  PROBLEM 


The  model  program  was  received  from  NOAA  in  printout 
form  and  was  typed  onto  my  *A  disk  on  the  NPS  system  with 
each  subroutine  filed  as  a  separate  Fortran  program.  During 
this  phase,  it  became  evident  that  the  full  60  bit  capa¬ 
bility  of  the  CDC  computer  had  been  utilized  .  To  achieve  a 
comparable  level  of  accuracy  on  the  IBM  system,  it  was 
necessary  to  use  double  precision.  Rather  than  change  all 
the  variables  to  double  precision  before  more  familiarity 
was  gained  with  the  program,  it  was  decided  to  use  the 
autodbl  function  which  would  automatically  change  all  the 
variables  and  functions  to  IBM  double  precision.  Mr. 
Hilleary  at  the  computer  center  was  most  helpful  in  this 
endeavor. 

Because  of  the  complexity  of  the  program  and  its  modular 
structure,  it  was  decided  that  a  ground-up  module  test 
method  would  prove  most  advantageous.  Thus  once  all  the 
subroutines  were  copied  and  checked,  a  short  HATFIV  program 
was  written  around  each  subroutine  to  test  its  correct  oper¬ 
ation.  The  WATFIV  compiler  was  utilized  because  of  its 
excellent  error  messages.  Each  program  supplied  the 
required  input  variables,  called  the  respective  subroutine 
and  displayed  the  desired  variables  which  were  spot-checked 
using  a  handheld  calculator.  During  the  initial  parts  of 
this  phase,  r.o  attempt  was  made  to  integrate  the  subrou¬ 
tines,  i.e.  each  program  tested  one  subroutine  or  function 
only,  in  an  attempt  to  minimize  the  amount  of  debugging 
required.  However,  once  all  the  subroutines  had  been  tested 
individually,  some  integrated  testing  was  conducted.  The 
modular  testing  procedure  additionally  served  as  an  excel¬ 
lent  means  of  gaining  familiarity  with  each  part  of  the 


program.  This  proved  to  be  of  great  benefit  in  later  stages 
of  debugging  and  modification. 

Simultaneous  with  the  testing  and  debugging  process,  it 
was  realized  that  actual  rawinsonde  data  would  be  required 
once  the  program  was  ready  to  run.  To  avoid  a  lag  in  the 
development,  an  effort  was  made  to  procure  some  representa¬ 
tive  data  fairly  early.  The  dates  and  locations  were  chosen 
to  coincide  with  a  detailed  experiment  run  by  the  University 
of  Arizona  and  the  Air  Force  £Bef.  4]  at  Sunset  Canyon,  Az. 
and  Truckee,  Nev.  in  1977-78.  The  Naval  Environmental 
Prediction  and  Research  Facility  (NEPRF) ,  located  in 
Monterey,  Ca.,  provided  rawinsonde  data  upon  letter  request. 
The  rawinsonde  data  output  was  in  the  form  of  raw  tabulated 
values  for  the  mandatory  and  significant  reporting  levels 
and  a  summary  table  which  listed  the  desired  variables 
including  atmospheric  pressure,  temperature,  humidity,  wind 
speed  and  direction  and  absolute  altitude.  Additionally 
graphs  were  also  received  but  were  not  utilized.  The  summa¬ 
rized  values  were  filed  in  the  proper  format  in  the  DATA  I'! 
file. 

The  next  step  was  to  write  a  data  processing  subroutine 
that  would:  1.  process  the  raw  data  into  the  units  required 
by  the  subroutines  and,  2.  calculate  the  necessary  derived 
variables  such  as  the  humidity  gradient  and  the  wind  shear. 

Nith  the  completion  of  the  data  processing  subroutine, 
all  was  in  readiness  for  the  development  of  the  main 
program.  The  function  of  the  main  program  was  to  step  the 
calculations  from  the  first  data  input  set  (corresponding  to 
the  lowest  slab  of  a  given  atmosphere)  sequentially  through 
each  slab  until  the  last  data  set  was  reached.  The  main 
program  also  initialized  a  number  of  common  variables  that 
were  used  in  various  subroutines  throughout  the  program, 
lastly,  it  served  as  a  framework  for  the  proper  sequencing 
of  the  subroutine  calls. 


As  might  be  expected  in  complex  programs,  even  with  the 
extensive  modular  testing  previously  conducted,  a  signifi¬ 
cant  amount  of  debugging  was  required  once  the  program  was 
finally  integrated.  When  attempts  to  decipher  the  Fortran 
error  messages  proved  fruitless,  the  WATFIV  compiler  was 
again  utilized.  The  main  difficulty  in  debugging  was  trying 
to  trace  the  errors  through  a  series  of  imbedded  subrou¬ 
tines.  For  example  one  of  the  most  difficult  problems  to 
detect  at  this  point  was  caused  by  an  underflow  error.  The 
EXPNX  function  was  written  to  prevent  underflow  or  overflow 
errors  from  occurring.  As  received,  it  tested  to  see  if 
variables  exceeded  the  limits  of  the  CDC  computer  (approxi¬ 
mately  1.E±300) .  During  the  early  testing,  it  was  realized 
that  these  limits  would  have  to  be  changed  to  correspond  to 
IBM  system  3033*s  capacity  (approximately  1.5±75)  .  If  the 
variable  in  question  exceeded  these  limits,  the  function 
EXPNX  would  revalue  the  variable  to  the  respective  limit 
imposed  by  the  computer.  What  was  not  immediately  apparent 
however,  was  that  the  corrected  variable  was  subsequently 
multiplied  by  another  small  number  (on  the  order  of  1.  5-3) 
in  another  subroutine  which  ultimately  caused  the  underflow. 
While  in  itself,  the  error  was  relatively  easy  to  detect,  it 
proved  difficult  and  time  consuming  to  trace  since  the 
offending  subroutine  was  buried  in  other  subroutines. 

Cne  technique  that  also  proved  effective  was  to  draw  a 
large  flowchart  which  listed  each  subroutine  and  function 
contained  in  the  program.  Additionally  each  subroutine  or 
function  was  annotated  with  the  other  subroutines  that 
called  it  as  well  as  those  that  were  called  by  it.  This  was 
very  useful  in:  1.  gaining  the  ’big  picture’  of  the  struc¬ 
ture  of  the  program,  2.  checking  that  the  required  vari¬ 
ables  were  indeed  entered  or  printed  and  that  each  was 
listed  in  the  call  statement,  3.  detecting  several  vari¬ 
ables  that  had  inadvertently  undergone  name  changes  in  the 


genesis  of  the  prograa  at  NOAA  and  had  subsequently  becoae 
sources  of  confusion  in  the  translation,  and  4.  tracing 
suspect  variables  back  to  their  origin. 


III.  SOLUTION 


A.  THEOBI 


VanZandt  et  al,  [Bef.  3],  assume  the  existence  of  thin 
horizontal  layers  of  turbulence  in  an  otherwise  non- 
turbulent  atmosphere.  These  layers  are  evident  in  radar 
[Bef.  3]  and  direct  balloon  measured  Cn2  profiles,  [Bef.  5] 
and  [Bef.  6].  They  vary  with  time  and  location,  though  they 
are  generally  of  limited  horizontal  area.  The  large  scale 
turbulent  eddies  in  this  thin  region  are  inherently  inhomo¬ 
geneous,  anisotropic  and  non-steady.  Nevertheless,  within 
the  turbulent  layers,  the  turbulence  approaches  homogeneous 
and  isotropic  conditions  if  the  scale  length  is  taken  small 
enough  (i.e.  scale  lengths  much  shorter  than  the  thickness 
of  the  layer)  .  Furthermore,  it  is  assumed  that  all  the 
different  sized  eddies  come  into  existence  within  a  short 
time  span,  relative  to  the  life  span  of  the  turbulent  eddy 
itself.  Hith  these  approximations,  the  velocity  fluctua¬ 
tions  can  be  considered  to  be  steady-state  and  the  turbu¬ 
lence  structure  parameter,  Cn2,  can  be  estimated  for  the 
turbulent  layer.  To  estimate  the  turbulence  throughout  an 
entire  vertical  shaft  of  atmosphere,  it  is  also  necessary  to 
estimate  the  fraction,  F,  of  the  entire  shaft  that  is  turbu¬ 


lent  , 


The  calculation  and  combination  of  these  two  esti¬ 


mates  form  the  basis  for  this  model  and  the  program. 

For  the  case  of  homogeneous,  isotropic  turbulence, 
Tatarski  [Ref.  1],  estimates  Cn2  by; 


Cn2  =  a2  •  ot*  •  l.  «/3  •  fjz 


(3.1) 
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where  a  is  a  universal  constant  generally  taken  to  be  2.8, 
«**  is  the  ratio  between  the  eddy  viscosity  and  the  eddy 
thermal  diffusivity  generally  taken  to  be  1.3,  L#  is  the 
outer  scale  length  of  the  turbulence  and  M  is  the  vertical 
gradient  of  the  refractive  index. 

lo  is  assumed  to  be  of  the  order  of  the  thickness  of  the 
turbulent  layer  itself.  An  expression  for  L« ,  also  given  by 
Tatarski,  is; 

L.  =  {  e  /  S3  )  V2  (3.2) 

where  e  is  the  dissipation  rate  of  turbulent  kinetic  energy 
and  S  is  the  vertical  gradient  in  the  horizontal  win  I 
velocity  (i.e.  J  V/  J  z)  .  Bather  than  try  to  measure  these 
values,  the  developers  of  this  model  selected  a  value  of  Lo 
that  would  normalize  the  theoretical  Cn2  to  that  produced  by 
their  experimental  radar  measurements.  This  value  of  L. 
turned  out  to  be  10  meters. 

The  expression  for  M  is  given  by  [Hef.  3]  as; 

M  =  -77.6  •  I.E”6  •  (  p  •  dine  /  liz  )  (3.3) 

•  (  1  ♦  15500  g/T  •  [  1-.  5  (  &ln  g/a  z)  /  (  din  6  /  dz)  ] 

where  p  is  the  atmospheric  pressure  in  millibars,  T  is  the 
absolute  temperature  in  Kelvin,  6  is  the  potential  tempera¬ 
ture  in  Kelvin,  g  is  the  specific  humidity  and  z  is  the 
altitude  in  meters.  To  calculate  fl,  the  vertical  resolution 
of  the  meteorological  data  should  be  on  the  order  of  the 
layer  thickness,  Lo  .  Current  rawinsonde  measurements  are 
somewhat  coarser,  on  the  order  of  100  meters  (for  routine 
equipment,  special  equipment  is  available  which  is  capable 
of  producing  a  much  finer  dataset)  ,  however  the  rawinsonde 
data  is  frequently  the  best  data  available.  Using  egns. 
(3.3)  and  (3.1),  it  is  possible  to  estimate  the  value  of  Cn2 
in  the  turbulent  layers  from  routine  rawinsonde  data. 


A  slab  of  atmosphere  on  the  order  of  10  meters  thick 
will  usually  contain  at  least  one  set  of  turbulent  and  non- 
turbulent  layers.  Having  previously  estimated  the  value  of 
Cn2  in  the  turbulent  layer#  the  average  Cn2  for  slab  i  is; 


Cn2  (i)  =  Cn2  (i)  •  F  (i)  (3.4) 

where  F(i)  is  the  mean  fraction  of  slab  (i)  that  is  turbu¬ 
lent.  There  are  no  expressions  to  solve  for  F  directly  so  a 
statistical  distribution  is  used. 

First  one  calculates  the  wind  shear,  which  is  the 
difference  between  the  horizontal  wind  velocity  vector  at 
two  consecutive  altitudes  (note  that  any  vertical  wind 
component,  such  as  due  to  convective  flow,  is  neglected)  . 
At  any  point,  the  total  shear,  St,  can  be  thought  of  as 
being  composed  of  two  parts,  the  mean  shear,  Sm,  and  a  fluc¬ 
tuating  part,  Sf.  Then  at  any  point; 

St  =  Sm  +  Sf  (3.5) 

The  coarseness  of  the  rawinsonae  measurements  will  provide 
sufficient  accuracy  for  the  calculation  of  Sm,  but  will 
preclude  the  direct  calculation  of  Sf  which  would  require 
the  separation  of  data  points  to  be  on  the  order  of  several 
meters.  Thus  it  was  assumed  that  Sf  has  a  Gaussian  distri¬ 
bution  Z  (  )  and  a  standard  deviation,  sigma,  which  is  inde¬ 
pendent  of  Sm  so  that; 

Z  (Sf/sigma)  =  exp  (  -Sf  2/2«sigma2)  /  (2«pi)1/2  (3.6) 

Then  the  probability,  f,  of  a  total  shear  in  the  increment 
St  to  St  +  dSt  is; 

f(St)  dSt  =  Z[  (St-Sm)  /sigma  ]  d(St/sigma)  (3.7) 
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It  was  also  assumed  that  a  region  is  turbulent  only  if  St 
exceeds  a  critical  shear.  Sc.  Then  F  is  given  by; 

|SC| 


F  ( |  taucrit] ,  ]  taumean] ) 


If  (St)dSt  +  / f  (St)dSt 

>|SJ 

=  P  {-! taucrit j - Jtaumean | )  +  Q ( | taucrit |- 1 taumean! ) 


k 


(3.8) 


where  |  taucrit  |  =  |Sc|/sigma,  |  taumean  |  =  |Sm| /sigma,  P(  ) 
is  the  cumulative  distribution  function  of  Z (  )  and  Q{  )  =  1 
-  P(  ).  Thus,  F  varies  with  1 taumean! ,  the  magnitude  of  the 
wind  shear  and  inversely  with  | taucrit |,  the  atmospheric 
stability. 

From  eqn.  (3.8) ,  F  is  a  function  of  two  variables, 
taucrit  and  taumean,  which  in  turn  are  functions  of  three 
variables.  Sc,  Sm  and  sigma.  Sm,  as  mentioned,  is  computed 
from  the  rawinsonde  data.  Sc  is  calculated  from; 

Sc  =  [  (g  •  dlnfl/}z)  /  He  (3.9) 

=  2  (  g  •  fcln  8  /  dz  ) 

The  critical  value  for  Richardson's  number,  Rc,  for  the 
onset  of  turbulence  is  taken  to  be  .25  and  this  value  was 
substituted  in  the  second  equality  of  eqn.  (3.9).  The  value 
of  sigma  is  taken  as  .01/s  in  the  troposphere  and  .015/s  in 
the  stratosphere,  [Bef.  3].  Fith  these  variables  in  hand, 
F  (i)  for  each  slab  can  be  calculated  from  eqn.  (3.8)  and 
ultimately  Cn2  from  eqn.  (3.4). 

B.  HETHOD 

The  program  as  received  was  similar  to  IBM  Fortran  IV 
overall,  although  there  were  some  significant  differences. 
The  CDC  to  IBM  Conversion  Guide,  [Ref.  7],  proved  most 
helpful  as  did  the  staff  at  the  Computer  Center  and  local 
expertise  in  the  Aeronautical  Engineering  Dept,  at  NPS.  The 


most  frequent  modifications  that  were  made  have  been  listed 
in  Table  1  . 


TABLE  1 

CDC  to  IBM  Translation 


Item 

CDC  code 

IBM  translation 

Comment  card 

uses  *  or  C 

substitute  C  for  * 
in  col.  #1 

Continuation 

card 

uses  *  1*  in 
ccl.  #6 

substitute  *  for  1 

Statement 

separation 

uses  A=1$B=2$C=3 
for  example 

write  each  statement 
on  different  lines, 
i.e.  A=1 

3=2,  etc. 

Initiali ze 
statements 

uses  A=B=C=0 

write  each  statement 
or  equality  on 
different  lines, 
i.e.  A  =0 

B=0,  etc. 

Program 

statements 

uses  Program 
statement  at  the 
start  of  the  main 
subroutine 

delete  Program 
sta tement 

Name  length 
of  variable 

permits  up  to  7 
characters  in 
name 

rename  variables  so 
that  maximum  length 
doesn’t  exceed  6 
characters 

Name  length 
of  subroutine 

same  as  above 

same  as  above 

Delimiters 

uses  PRINT,,ABC" 
for  example 

substitute 
*  in  place  of  " 

Print 

statements 

uses  PEI  NT*  A, B, C 
or  PRINT  100, X 

substit  ute 

WRITE  (5  ,  100)  X 
as  required 

Integer 

division 

uses  LAST=N QOAD/2 
where  LAST  and 

HQUAD  are  integers 
and  NQOAD  is  odd 

rewrite  using 
different  approach 

Once  the  program  was  running  satisfactorily,  the  Disspla 
graphics  routine  was  developed.  Since  the  program  was 
written  as  single  precision  and  then  converted  into  double 


precision  through  the  use  of  the  autodbl  function,  all  of 
the  output  values  were  necessarily  of  double  precision  type 
and  it  was  impossible  to  selectively  override  the  function 
so  as  to  declare  some  of  the  variables  as  single  precision. 
Because  Disspla  will  not  accept  double  precision,  two  alter¬ 
natives  presented  themselves  at  this  point,  either;  1. 
write  a  separate,  single  precision  plotting  program  that 
would  read  the  main  program’s  data  output  file  (hereafter 
the  DATAOOT  file),  as  single  precision,  and  plot  the  values 
or,  2.  convert  the  entire  program  to  double  precision  except 
for  two  single  precision  variables  which  would  serve  as  the 
x  and  y  coordinates  for  the  Disspla  routine.  Since  the 
program  for  alternative  #1  was  already  working,  it  seemed 
the  most  expedient  choice  and  therefore  a  separate,  short 
plotting  program  was  developed.  Once  it  was  running 
however,  it  became  apparent  that  having  to  run  first  one 
program  and  then  another  was  overly  cumbersome,  especially 
for  a  program  to  be  made  available  for  general  use. 
Alternative  #2  was  then  developed,  i.  e.  all  the  program 
variables  were  declared  as  REAL*8  except  for  two  single 
precision,  dummy  variables,  PLOTAX  and  PLOTAY,  and  the  plot¬ 
ting  routine  was  incorporated  as  a  section  of  the  main 
program.  This  proved  very  effective  because  the  user  now 
received  the  graphics  at  a  Tektronics  618  terminal  followed 
by  the  tabulated  data  in  a  DATAOUT  file  which  is  available 
for  browsing  or  printing  from  the  user’s  Filelist  (FLIST)  . 

C.  INITIAL  CALCULATIONS 

The  main  program  first  declares  all  the  variables  and 
initializes  the  necessary  constants  which  are  used  in  the 
program  and  reads  the  program  options  (i.e.  whether  or  not  a 
plot  is  desired  as  well  as  the  number  of  runs  desired) . 
Next  the  program  reads  the  number  of  data  sets,  N  (i.e.  the 
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number  of  levels  reported  by  the  rawinsonde) ,  from  the 
DATAIN  file.  The  variable  format  is  15.  Following  that  a 
'Do  loop',  iterated  N  times,  reads  the  data  sets  from  the 
DATAIN  file.  For  each  level,  the  data  set  format  is 
5(F10.5),  F10.2  and  the  variables  are  in  the  following 

order:  pressure  (millibars),  temperature  (Celsius),  rela¬ 

tive  humidity  (55/100),  wind  speed  (meters/  second),  wind 
direction  (compass  heading  in  degrees  relative  to  true 
North)  and  absolute  altitude  (meters).  The  respective  vari¬ 
able  names  are:  PRESS,  TEMP,  HOMID,  RS,  RD,  and  Z. 

The  data  processing  subroutine  (DATPRO)  is  called  with 
these  six  variables  plus  N.  The  subroutine  calculates; 

TEMPA  (the  absolute  temperature,  TEMP  +  273.) 

T2KARY  (the  thickness  of  the  slab,  Z(I+1)-Z(I)) 

DQDZRY  (the  average  gradient  of  the  humidity  through  the 

slab) 

PARY  (the  weighted  average  of  the  pressure  for  the  slab) 

TAEY  (the  weighted  average  of  the  temperature  for  the 

slab) 

CARY  (the  weighted  average  of  the  humidity  for  the  slat) 

PARY,  TARY  and  QARY  were  originally  taken  as  simple  aver¬ 
ages,  but  were  changed  to  weighted  averages  in  order  to 
smooth  the  roughness  introduced  when  point  rawinsonde  values 
are  used  for  continuously  varying  (spatially  and  temporally) 
physical  parameters.  The  generic  weighting  equation  used 
was ; 

X  (I)  =  (  Y  (I - 1 )  +  3.  •  Y  (I)  (3.10) 

♦  3.  •  Y  (1  +  1)  +  Y  (1+2)  ]  /  8. 
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The  value  of  THETA,  the  potential  temperature,  is  calcu¬ 
lated  next  from; 

THETA  (I)  =[  (PRESS  (1) /PRESS  (I)  )**  (R/Cp)  ]«TEMPA  (3.11) 

where  R  is  the  gas  constant  and  Cp  is  the  specific  heat. 

The  wind  vector  is  then  broken  down  into  the  usual  u  and 
v  components  (OX  and  OX  respectively)  and  the  wind  shear  is 
calculated  from; 

SHEAR  (I)  =  [  (OX  (1  +  1)  -  OX  (I)  )2  +  (3.12) 

(  OX  (1+1)  -  OX  (I)  )2  31/2  /  THRARX  (I) 

Next,  the  stability  parameter,  STBARX,  is  calculated  by; 

STBAEX(I)  =  G«[  LOG  (THETA  (1+  1)  )  -LOG  (THETA  (I)  )  ]  (3.13) 

/  C  Z  (1+1) -Z  (I)  ] 

linally  Richardson’s  number,  RIARX,  for  the  slab  is 
calculated  by; 

RIARX  (I)  =  STB  AR  X  ( I)  /  (SHEARX(I))2  (3.14) 

The  values  of  TARX,  PARX,  QARX,  STBARX,  SHEARY,  DQDZRY, 
RIARX  and  THRARX  are  all  arrays  with  N-1  elements.  These 
arrays  are  returned  to  the  main  program.  The  primary  ’Do 
loop’  of  the  main  program,  which  iterates  over  the  altitudes 
corresponding  to  the  input  data,  then  calls  the  subroutines 
SHRSGM,  STESGM  and  INTLGS .  INTLGS  calls  each  of  the 
remaining  subroutines  and  functions.  A  flowchart.  Figure 
3.1,  depicts  the  general  flow  within  the  program.  Rather 
than  describe  each  subroutine,  a  brief  summary  of  each  is 
included  as  Appendix  B. 
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The  end  product  is  the  value  of  PCNSQ  (or  Cn2)  calcu¬ 
lated  for  each  slab.  The  average  altitude  of  the  slab, 
MIDA1T,  is  used  for  the  graph  and  tabulation.  The  Cn2 
values  with  the  corresponding  MIDALT's  are  stored  in  two  N-1 
element,  single  precision  arrays,  PLOTAX  and  PLOTAY  respec¬ 
tively.  At  this  point,  if  the  number  of  runs  is  specified 
as  '2'  in  the  DATAIN  file,  the  program  reads  in  all  the 
datasets  again  and  recalculates  the  Cn2  profile  (this  option 
will  be  further  discussed  in  Chap.  4) .  Next,  if  specified 
(i.e.  if  PLOT  =  '2*  in  the  DATAIN  file),  the  Disspla  routine 
is  invoked  which  plots  the  value  of  log  Cn2  on  the  x-axis 
vs.  HIDAIT  on  the  y-axis.  The  graph  displayed  is  available 
in  hardcopy  form  at  the  Tektronics  613  terminal.  Finally, 
the  Cn2  and  MIDALT  values  are  also  listed  in  the  TURB  2 
DATAOUT  file  which  has  been  transferred  to  the  user's  'A' 
disk.  This  file  may  then  be  browsed  or  printed  as  desired. 

D.  I1P0T  FOSHAT 

The  following  section  describes  the  step-by-step  proce¬ 
dure  to  be  followed  when  executing  TURB2.  Prior  to  execu¬ 
tion  of  course,  it  is  necessary  to  obtain  copies  of  the 
required  programs,  namely  TURB  DATAIN,  T0R32  FORTRAN  and 
optionally,  REFP10I  CATAIN.  Additionally,  an  exec,  TURB 2 
EXEC,  is  available  and  is  very  convenient  to  use. 

The  first  step  is  to  type  the  rawinsonde  data  into  the 
TURB  DATAIN  file.  The  easiest  method  is  to  copy,  rename  and 
write-over  an  old  DATAIN  file  (see  Appendix  A),  however  the 
format  is  listed  in  Table  2  . 

When  the  rawinsonde  data  is  received,  each  data  set 
(i.e.  a  set  of  all  the  measurable  variables  at  a  given  alti¬ 
tude)  will  fall  into  one  of  three  categories,  either; 

1.  the  dataset  has  balloon  measured  data  (pressure, 
temperature  and  humidity  along  with  radar  altitude)  tut 


TABLE  2 

TOBB  DATA IN  File  Format 


Card  (s) 

Cols. 

Type 

Variable 

Remarks 

1 

20 

15 

PLOT 

enter  ’ 1  *  if  a  plot  is 
desired,  *2’  if  not 

2 

20 

15 

#RUNS 

enter  the  number  of 
runs  desired  (1  or  2) 

3 

16-20 

15 

N 

number  of  datasets  to 
be  entered 

4->N+4 

1-10 

F10.5 

PRESS 

atmospheric  pressure 
(mi Hi  bars) 

ft 

11-20 

F10.5 

TEMP 

temperature  (degrees 
Celsius,  not  absolute) 

It 

21-30 

F10.5 

HUMID 

relative  humidity  (?) 

ft 

31-40 

F10.5 

WS 

wind  speed  (meters/sec) 

ft 

41-50 

F10.  5 

WD 

wind  direction  (compass 
direction,  degrees) 

ft 

51-60 

F10.2 

Z 

absolute  altitude 

(meters) 

Note:  If  #BUNS  =  1.  the  DATAIN  file  is  complete.  If  ftRUUS  =  2, 
copy  cards  3  thru  N+4  and  insert  them  immediately  following 
card  N+4  (see  Appendix  A)  . 

is  missing  the  ground  radar  generated  data  (wind 
velocity)  , 

2.  the  radar  wind  speed  and  direction  data  appears  but 
the  balloon  data  is  missing,  or 

3.  all  data  is  recorded  (at  the  7  or  8  mandatory  levels 
only)  . 

Typically,  either  balloon  or  radar  data  is  taken  at  80  to 
100  points. 

Note  that  the  Datain  file  and  the  program  accept  only 
complete  data  sets.  Enter  a  value  of  ’999.9*  for  any  data 
points  that  are  missing.  The  program  will  then  do  a  semi¬ 
log  interpolation  for  the  missing  values  using  one  of  two 
methods  as  follows; 


1.  if  IRONS  =  1,  all  the  missing  values  are  assigned 
interpolated  values.  In  this  case,  the  program  calcu¬ 
lates,  plots  and  tabulates  only  one  Cn2  profile  or 

2.  if  #RUNS  =  2,  the  program  rereads  the  raw  datasets  and 
interpolates  only  for  missing  values  of  WS  and  WD.  Thus 
those  datasets  missing  PR SSS,  TEMP  and  HUMID  are  deleted. 
In  this  case,  a  new  Cn2  profile  is  estimated,  which  is 
plotted  and  tabulated  in  addition  to  the  previous  one 
(i.e.  when  IfiUNS  =  1).  Setting  IRONS  -  2  nearly  doubles 
the  number  of  executions  and  thus  the  run  time.  The  use 
of  this  option  will  be  discussed  in  Chap.  4. 

The  program  may  now  be  compiled  (if  necessary)  and 
executed,  it  does  not  require  batch  processing  and  the 
output  should  be  ready  within  several  minutes.  Note  that  it 
is  possible  to  store  many  data  sets  in  one's  file,  but  the 
TURB2  EXEC  reads  from  TURB  DATAIN.  Thus  the  user  may  find 
it  necessary  to  either  rename  his/her  DATAIN  file  or  change 
the  EXEC. 

If  the  user  wishes  to  display  a  second  Cn2  plot  simulta¬ 
neously  with  the  model's  estimated  plot,  the  values  for  the 
second  plot  are  stored  in  a  second  DATAIN  file  called 
EEFPICT  DATAIN.  This  is  not  the  same  as  setting  IF.UNS  =  2. 
In  this  case  a  previously  measured  or  estimated  Cn2  profile 
is  input  and  the  program  merely  plots  the  input  values.  As 
in  the  previous  case,  the  easiest  method  is  to  copy  and 
rewrite  an  old  file,  however  the  format  is  as  follows: 


TABLE  3 

Reference  Plot  Da tain  File  Format 


Card  (s) 


Cols. 


Type 


Variable 


Remarks 


If.  BESPITS 


1.  C  OH  PI  B I  SOI 

Once  the  program  was  debugged  and  executing  reliably,  it 
was  necessary  to  compare  the  output  to  that  produced  by  the 
original  program  at  NOAA.  This  procedure  served  as  an 
acceptance  test  to  validate  the  operation  of  the  translated 
program.  To  this  end,  a  random  input  and  output  dataset  was 
requested  from  NOAA.  When  the  translated  program  was  run 
with  the  test  input  on  the  NPS  computer  system,  the  output 
was  found  to  agree  closely  with  output  produced  by  the  orig¬ 
inal  code  executed  on  NOAA's  CDC  computer  (see  Tables  4  and 
5,  Appendix  D)  .  Table  4  was  produced  using  the  first  method 
described  (i.a.  #RtJNS  =  1)  while  Table  5  was  produced  using 
the  second. 

As  can  be  seen  from  Figure  4.1  and  Table  4,  when  the 
first  interpolation  method  was  used  (i.e.  interpolated 
values  are  substituted  for  all  missing  datapoints) ,  there 
are  noticeable  disagreements,  but  overall  the  I3M  output  is 
within  17.93b  of  the  NOAA  results  (if  one  neglects  the  large 
spikes)  .  The  difference  in  the  accuracy  between  low  and 
high  altitudes  was  insignificant.  There  is  one  large  excur¬ 
sion  around  15  km.  that  is  without  explanation.  Also,  there 
were  several  spikes  in  the  NOAA  data  at  2  km.  and  at  6.3  km. 
that  seem  suspect  and  were  not  repeated  in  the  IBM  output. 

As  can  be  seen  from  Figure  4.2  and  Table  5,  when  the 
second  interpolation  method  was  used  (i.e.  interpolated 
values  are  substituted  for  datasets  missing  RS  and  WD  only, 
with  the  other  datasets  being  deleted) ,  the  output  tended  to 
be  in  closer  agreement  with  the  NOAA  results.  However,  in 
several  areas,  the  IBM  data  tended  to  respond  more  slowly 
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Figure  <1.1  Comparison  of  Results,  IRONS 


3anxinv  sa  s..no  ooi 


Comparison  of  Besults,  #BUNS 


than  the  NOAA  results  (because  of  the  fewer  datasets)  and 
some  of  the  spikes  in  the  NOAA  data  were  ’smoothed*  into 
insignificance.  The  underlying  question  is  whether  the 
resolution  of  the  rawinsonde  data  is  sufficient  to  predict 
the  variations  in  the  value  of  Cn2  at  the  scale  shown  or 
conversely,  is  the  data  being  overprocessed  to  the  point 
that  the  spikes  are  merely  ’noise’.  To  answer  the  question, 
it  will  be  necessary  to  process  more  datasets  and  also  to 
measure  the  Cn2  profile  experimentally  so  that  comparisons 
can  be  made. 

Overall,  with  #RUNS  =  2,  the  IBM  output  was  found  to  be 
within  14.75S  of  the  NOAA  results.  In  the  region  below  12 
km.,  the  agreement  was  within  20.1%  of  the  reference  data. 
The  spike  at  6575  meters  seems  to  be  produced  by  a  glitch  in 
the  program  (because  the  wind  shear  is  zero  for  several 
consecutive  slabs)  and  probably  doesn’t  predict  the  actual 
Cn2  accurately.  Above  12  km.,  the  agreement  is  excellent, 
averaging  5.0%  and  without  any  significant  excursions. 

There  are  two  primary  reasons  why  the  output  doesn’t 
agree  exactly,  as  one  might  expect.  The  first  is  due  to  the 
fact  that  the  data  processing  subroutine  used  in  the  NOAA 
program  was  not  received  with  the  rest  of  the  code  and  thus 
a  substitute  had  to  be  written.  The  original  (i.e.  NOAA) 
data  processing  routine  was  received  later  with  the  compar¬ 
ison  test  data.  Although  the  two  are  very  similar,  there 
are  some  differences  and  thus  slightly  different  values  are 
input  to  the  main  program  in  each  case.  The  second  reason 
is  that  the  NOAA  program  has  undergone  some  modifications  in 
the  time  between  the  receipt  of  the  original  code  (from 
which  the  IBM  code  derived)  and  the  receipt  of  the  compar¬ 
ison  test  data.  While  it  is  highly  unlikely  that  any  large 
differences  would  be  introduced  by  the  modifications,  it  is 
almost  certain  that  some  small  differences  crept  in.  Since 
the  NOAA  program  itself  is  still  undergoing  development,  it 
is  impractical  to  try  to  keep  current  translated  copies. 
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However,  with  these  differences  in  mind,  the  output  of 
the  translated  code  is  quite  satisfactory  for  the  purposes 
specified  and  has  the  potential  for  being  much  superior  to 
the  accuracy  of  any  of  the  previously  available  models. 
Also  the  ease  of  operation  and  useful  graph/  tabular  output 
format  make  it  a  very  practical  and  useful  addition  to  the 
NPS  library. 
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7.  CONCLUSIONS  AND  BECO MMEHDATIOHS 

The  NO AA  program  differs  primarily  in  the  smoothing 
algorithm  used  that  produces  a  different  dataset  for  the 
main  program.  The  translated  program  was  modified  so  that  a 
comparison  could  be  made  between  the  two  methods  of  data 
processing  described  earlier.  The  intent  of  the  comparison 
was  to  see  if  the  first  method  was  actually  more  accurate  or 
if  the  larger  number  of  datasets  only  tended  to  increase  the 
•noise* . 

While  the  results  indicated  that  the  first  method  is 
better,  it  is  difficult  to  make  a  judgement  based  on  the 
results  of  a  single  dataset.  However,  for  routine  use,  the 
first  method  is  recommended  (i.e.  set  #BUNS  =  1  in  the 
DATAIN  file)  . 

One  of  the  strengths  of  this  model  is  that  it  requires 
only  readily  available  rawinsonde  data  to  be  input,  but  that 
simplicity  is  also  a  basic  source  of  error  that  may  occur  in 
three  forms. 

The  first  error  comes  about  because  the  rawinsonde  data, 
while  available  for  a  great  many  locations  and  times,  is  not 
available  for  all  places  and  times.  If  one  wishes  to 
compare  the  model  predictions  with  measurements  from  a  radar 
or  scintillomete r,  it  is  necessary  to  assume  that  the  atmos¬ 
pheric  parameters  recorded  by  the  rawinsonde  are  approxi¬ 
mately  eguivalent  to  those  present  at  the  measurement  site. 
This  may  or  may  not  be  a  reasonable  assumption  because  of 
the  spatial  or  temporal  variation  in  the  parameters.  Tor 
short  distances  (less  100  miles  and  without  large  variation 
in  geographical  features  like  mountains  or  oceans) ,  and 
times  (2-4  hours  depending  on  the  time  of  day) ,  the  assump¬ 
tion  is  probably  reasonable.  Obviously,  as  distance  or  time 


less  correlation 


from  the  observation  point  increases, 
between  the  model  and  the  measured  data  are  expected. 

The  second  error  derives  from  the  fact  that  scintillom¬ 
eter  measurements  are  time  averaged  line  integrals,  whereas 
the  raw  insonde  balloon  carries  an  array  of  point  sensors. 
Also  the  balloon  rises  over  a  one  to  two  hour  period  along 
an  oblique  path  as  a  function  of  the  wind  velocity,  while 
the  experimental  measurements  are  usually  made  along  a 
vertical  one.  For  these  last  two  reasons,  even  if  the 
balloon  launch  site  and  the  radar  or  scintillometer  were 
co-located  and  readings  were  made  simultaneously,  there 
would  still  be  some  residual  data  errors. 

The  third  source  of  error  occurs  only  if  there  are 
convective  layers  present  in  the  atmosphere.  Because  of  the 
difficulty  in  modeling  the  effects  of  such  layers,  the 
present  model  neglects  them. 

tfith  these  limitations  in  mind,  the  model  is  still  the 
best  available  and  should  be  of  utility  to  those  whose 
research  requires  the  estimation  of  Cn2(z).  As  the  parent 
program  continues  to  be  improved,  it  may  be  worthwhile  at 
some  future  time  to  add  the  corresponding  updates  to  this 
version  also.  Additionally,  I  think  that  the  program  could 
be  integrated  into  a  propagation  model,  such  as  G.U.T. S. , 
for  a  significant  increase  in  accuracy  over  the  present 
model. 


APPENDII  A 

TDHB  DAI  A  IN  FILE  FOBHAT  EXAMPLE 

This  appendix  contains  the  dataset  used  for  the  compar¬ 
ison  test.  It  is  also  to  be  used  as  an  example  of  the 
format  required  by  the  main  program.  Note  that  the  file 


contains 

two  copies 

of  the  same 

dataset 

.  This  is 

the  format 

used  when 

#R0NS  = 

2.  If  only  one 

run  is  desired,  the 

repeated 

dataset  may 

de  deleted  (i.e. 

there  would  only  be 

118  lines 

vice  236) . 

PLOT?  1 

(1-yes, 

2-no) 

#RUNS  =  2 

(1  or  2 

only) 

#  DATA  SETS  =  118 

PRESS. 

TEMP. 

HUMID. 

KS 

WD 

ALTITUDE 

(mbars) 

{  C) 

(S/100) 

(a/s) 

(degrees) 

(meters) 

834.9 

22.4 

.291 

4.1 

120. 

16  10. 

828.0 

20.2 

.261 

999.9 

999.9 

1680. 

999.9 

999.9 

999.9 

5.1 

158. 

1930. 

778.0 

15.5 

.357 

999.9 

999.9 

2200. 

999.9 

999.9 

999.9 

4.6 

164. 

2240. 

750.0 

12.0 

.424 

999.  9 

999.9 

2520. 

999.9 

999.9 

999.9 

5.  1 

173. 

2670. 

700.0 

6.9 

.575 

999.9 

999.9 

3100. 

999.9 

999.9 

999.9 

5.1 

176. 

3140. 

999.9 

999.9 

999.9 

5.  1 

186. 

3560. 

648.0 

.  8 

.883 

999.9 

999.9 

3720. 

999.9 

999.9 

999.9 

6.7 

212. 

3930. 

615.0 

-3.  1 

.963 

999.9 

999.9 

4140. 

999.9 

999.9 

999.9 

7.7 

243. 

4230. 

596.0 

-4.2 

.690 

999.9 

999.9 

4390. 

999.9 

999.9 

999.9 

7.7 

268. 

44  60. 
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579.0 

-6.  1 

.714 

999.9 

999.9 

999.9 

999.9 

999.9 

999.9 

554.0 

-8.5 

.664 

999.9 

999.9 

999.9 

529.0 

-10.1 

.372 

999.9 

999.9 

999.9 

513.0 

-12.  1 

.411 

999.9 

999.9 

999.9 

500.0 

-13.  8 

.619 

999.9 

999.9 

999.9 

483.0 

-16.0 

.532 

999.9 

999.9 

999.9 

460.0 

-17.2 

.054 

999.9 

999.9 

999.9 

452.0 

-18.8 

.304 

444.0 

-20.1 

.492 

999.9 

999.9 

999.9 

437.0 

-20.7 

.204 

999.9 

S99.9 

999.9 

410.0 

-24.5 

.247 

402.0 

-25.5 

.555 

400.0 

-25.7 

.417 

999.9 

999.9 

999.9 

389.0 

-27.0 

.041 

999.9 

999.9 

999.9 

999.9 

999.9 

999.9 

999.9 

995.9 

999.9 

340.0 

-35.  2 

.242 

999.9 

999.9 

999.9 

999.9 

99  9.9 

999.9 

316.0 

-39.7 

.314 

999.9 

999.9 

999.9 

300.0 

-42.0 

.000 

999.9 

999.9 

999.9 

999.  9 

999.9 

4620. 

7.7 

279. 

46  80. 

7.7 

278. 

4900. 

999.9 

999.9 

4960. 

7.2 

274. 

5140. 

999.9 

999.9 

5320. 

6.7 

276. 

5410. 

999.9 

999.9 

5550. 

6.7 

280. 

5690. 

999.9 

999.9 

57  50. 

6.7 

275. 

5980. 

999.9 

999.9 

6010. 

5.7 

266. 

6260. 

999.9 

999.9 

6380. 

4.  6 

256. 

6490. 

999.9 

999.9 

6510. 

999.9 

999.9 

6640. 

4.6 

256. 

6740. 

999.9 

999.9 

6760. 

5.  1 

257. 

6980. 

5.  1 

270. 

7230. 

999.9 

999.9 

7370. 

999.  9 

999.9 

74  10. 

4.  1 

278. 

7510. 

999.  9 

999.9 

7610. 

3.6 

266. 

7800. 

4.  6 

274. 

8120. 

4.6 

270. 

8430. 

999.9 

999.9 

8560. 

5.7 

259. 

87  20. 

6.2 

264. 

89  90. 

999. 9 

999.9 

9070. 

6.2 

271. 

9230. 

999.9 

999.9 

9420. 

8.  2 

280. 

9470. 
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296.0 

-42.6 

.000 

999.9 

999.9 

999.9 

283.0 

-41.9 

.000 

999.9 

999.9 

999.9 

999.9 

999.9 

999.9 

999.9 

999.9 

999.9 

250.0 

-46.0 

.000 

999.9 

999.9 

999.9 

237.0 

-48.5 

.000 

999.9 

999.9 

999.9 

999.9 

999.9 

999.9 

999.9 

999.9 

999.9 

999.9 

999.9 

999.9 

999.9 

999.9 

999.9 

200.0 

-53.5 

.000 

999.9 

999.9 

999.9 

185.0 

-55.  8 

.000 

999.9 

99  9.9 

999.9 

999.9 

99  9.9 

999.9 

999.9 

999.9 

999.9 

999.9 

999.9 

999.9 

150.0 

-56.3 

.000 

999.9 

999.9 

999.9 

999.9 

999.9 

999.9 

999.9 

999.9 

999.9 

130.0 

-5  9.6 

.000 

999.9 

999.9 

999.9 

999.9 

999.  9 

999.9 

999.9 

999.9 

999.9 

11^.0 

-58.0 

.000 

999.9 

999.9 

999.9 

108.0 

-59.4 

.000 

999.9 

999.9 

999.9 

100.0 

-58.3 

.000 

999.9 

999.  9 

999.9 

999.9 

999.9 

9510. 

11.3 

285. 

96  80. 

999.9 

999.9 

9810. 

14.4 

283. 

9930. 

21.6 

279. 

10220. 

23.2 

281. 

10500. 

999.9 

999.9 

10650. 

21.6 

282. 

10770. 

999.9 

999.9 

11000. 

24.  7 

278. 

11020. 

26.8 

273. 

112  80. 

27.3 

271. 

11540. 

26.8 

271. 

11790. 

27.3 

273. 

12050. 

999.9 

999.9 

12100. 

25. 2 

273. 

12320. 

22.7 

269. 

12600. 

25.7 

266. 

12880. 

27.8 

264. 

13160. 

27.3 

262. 

13430. 

25.2 

263. 

137  10. 

999.  9 

999.9 

13930. 

24.7 

265. 

13990. 

24.  2 

264. 

14280. 

25.7 

263. 

14570. 

999.9 

999.9 

14830. 

24. 7 

262. 

14860. 

23.  2 

264. 

15160. 

22.7 

264. 

15470. 

999.9 

999.9 

15500. 

20.6 

264. 

15730. 

19.0 

265. 

16000. 

17.0 

259. 

16300. 

999.9 

999.9 

16480. 

16.0 

266. 

16590. 
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999.9 

999.9 

999.9 

999.9 

999.9 

999.9 

999.9 

9  39.9 

999.9 

999.9 

999.  9 

999.9 

999.9 

999.9 

999.9 

999.9 

999.9 

999.9 

999.9 

999.9 

999.9 

70.0 

-6  1.0 

.000 

999.9 

999.9 

999.9 

999.9 

999.  9 

999.9 

999.9 

999.9 

999.9 

62.0 

-58.2 

.000 

60.0 

-58.5 

.000 

999.9 

999.9 

999.9 

999.9 

999.9 

999.9 

999.9 

999.9 

999.9 

999.9 

999.9 

999.9 

50.0 

-53.9 

.000 

999.9 

999.9 

999.9 

999.9 

999.9 

999.9 

999.9 

999.9 

999.9 

42.0 

-54.3 

.000 

999.9 

999.9 

999.9 

39.0 

-52.0 

.000 

999.9 

999.9 

999.9 

36.0 

-52.6 

.000 

999.9 

999.9 

999.9 

999.9 

99  9.9 

999.9 

999.9 

999.9 

999.9 

30.0 

-48.  2 

.000 

999.9 

999.9 

999.9 

28.0 

-47.0 

.000 

DATA  SE 

TS  =  118 

(delet 

14.4 

264. 

16870. 

13.9 

264. 

17150. 

13.4 

269. 

17430. 

12.4 

270. 

17710. 

10. 3 

266. 

17980. 

9.8 

266. 

18260. 

9.3 

266. 

18540. 

999.9 

999.9 

18710. 

8.2 

257. 

18820. 

8.2 

261. 

19100. 

5.7 

266. 

19380. 

999.9 

999.9 

19470. 

4.6 

262. 

19670. 

4.  1 

266. 

19960. 

4.6 

290. 

20240. 

6.7 

301. 

20520. 

4.6 

317. 

20800. 

999.9 

999.9 

20830. 

4.  1 

333. 

21120. 

2.1 

356. 

21440. 

1.0 

034. 

21760. 

999.9 

999.9 

21950. 

1.0 

064. 

22080. 

1 . 0 

267. 

22430. 

0.  5 

347. 

22770. 

999.9 

999.9 

22940. 

0.5 

007. 

231  10. 

0.5 

238. 

23450. 

0.  5 

229. 

23790. 

1.5 

229. 

24130. 

1.5 

207. 

24430. 

1.5 

195. 

24590. 

here  if 

IRONS  =  1) 

#  DATA  SETS 
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PRESS. 

TEMP. 

HO  HID. 

HS 

HD 

ALTITOD 

(abars) 

(  C) 

(*/100) 

(m/s) 

(degrees) 

(meters) 

834.9 

22.4 

.291 

4.1 

120. 

1610. 

828.0 

20.2 

.261 

999.  9 

999.9 

1680. 

999.9 

999.9 

999.9 

5.1 

158. 

1930. 

778.0 

15.5 

.357 

999.9 

999.9 

2200. 

999.9 

999.  9 

999.9 

4.6 

164. 

2240. 

750.0 

12.0 

.424 

999.9 

999.9 

2520. 

999.9 

999.9 

999.9 

5.1 

173. 

2670. 

700.0 

6.9 

.575 

999.9 

999.9 

3100. 

999.9 

999.9 

999.9 

5.1 

176. 

3140. 

999.9 

999.9 

999.9 

5.  1 

186. 

3560. 

648.0 

.  8 

.883 

999.9 

999.9 

3720. 

999.9 

999.9 

999.9 

6.7 

212. 

39  30. 

615.0 

-3.  1 

.963 

999.9 

999.9 

4140. 

999.9 

999.9 

^99. 9 

7.7 

243. 

4230. 

596.0 

-4.2 

.690 

999.9 

999.9 

4390. 

999.9 

999.9 

999.9 

7.7 

268. 

4460. 

579.0 

-6.1 

.714 

999.9 

999.9 

46  20. 

999.9 

999.9 

999.9 

7.  7 

279. 

4680. 

999.9 

999.9 

999.9 

7.7 

278. 

4900. 

554.0 

-8.5 

.664 

999.9 

999.9 

4960. 

999.9 

999.9 

999.9 

7.2 

274. 

5140. 

529.0 

-10.1 

.372 

999.9 

999.9 

5320. 

999.9 

999.9 

999.9 

6.7 

276. 

5410. 

513.0 

-12.  1 

.411 

999.  9 

999.9 

5550. 

999.9 

999.  9 

999.9 

6.7 

280. 

5690. 

500.0 

-13.8 

.619 

999.9 

999.9 

5750. 

999.9 

999.9 

999.9 

6.7 

275. 

59  80. 

483.0 

-16.0 

.532 

999.9 

999.9 

60  10. 

999.9 

999.9 

999.9 

5.7 

266. 

6260. 

460.0 

-17.2 

.054 

999. 9 

999.9 

6380. 

999.9 

999.9 

999.9 

4.6 

256. 

64  90. 

452.0 

-18.8 

.304 

999.9 

999.9 

6510. 

444.0 

-20.  1 

.492 

999.9 

999.9 

66  40. 
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999.9 

999.9 

999.9 

4.6 

256. 

6740. 

437.0 

« 

o 

<N 

1 

.204 

999.9 

999.9 

6760. 

999.9 

999.9 

999.9 

5.  1 

257. 

6980. 

410.0 

in 

• 

CM 

1 

.247 

5.  1 

270. 

7230. 

402.0 

-25.5 

.555 

999.9 

999.9 

7370. 

400.0 

-25.7 

.417 

999.9 

999.9 

74  10. 

999.9 

999.9 

999.9 

4.  1 

278. 

7510. 

389.0 

-27.0 

.041 

999.9 

999.9 

7610. 

999.9 

999.9 

999.9 

3.6 

266. 

7800. 

999.9 

999.9 

999.9 

4.6 

274. 

8120. 

999.9 

999.9 

999.9 

4.6 

270. 

84  30. 

340.0 

-35.2 

.242 

999.9 

999.9 

8560. 

999.9 

999.9 

999.9 

5.7 

259. 

8720. 

999.9 

999.9 

999.9 

6.2 

264. 

8990. 

316.0 

-39.7 

.  314 

999.9 

999.9 

9070. 

999.9 

999.9 

999.9 

6.2 

271. 

9230. 

300.0 

-42.0 

.000 

999.9 

999.9 

9420. 

999.9 

999.9 

999.9 

8.2 

280. 

9470. 

296.0 

-42.6 

.000 

999.  9 

999.9 

9510. 

999.9 

999.9 

999.9 

11.3 

285. 

96  80. 

283.0 

-41.9 

.000 

999.9 

999.9 

9810. 

999.9 

999.9 

999.9 

14.4 

283. 

99  30. 

999.9 

999.9 

999.9 

21.6 

279. 

10220. 

999.9 

999.9 

999.9 

23.2 

281. 

10500. 

250.0 

-46.0 

.000 

999.9 

999.9 

10650. 

999.9 

999.9 

999.9 

21.6 

282. 

10770. 

237.0 

-48.5 

.000 

999.9 

999.9 

11000. 

999.9 

999.9 

999.9 

24.7 

278. 

11020. 

999.9 

999.9 

999.9 

26.8 

273. 

11280. 

999.9 

999.  9 

999.9 

27.3 

271. 

11540. 

999.9 

999.9 

999.9 

26.8 

271. 

11790. 

999.9 

999.9 

999.9 

27.3 

273. 

12050. 

200.0 

-53.5 

.000 

999.9 

999.9 

12100. 

999.9 

999.9 

999.9 

25.2 

■  >. 

12320. 

185.0 

-55.8 

.000 

22.  7 

269. 

12600. 
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999.9 

999.9 

999.9 

25.7 

266. 

12880 

999.9 

999.9 

999.9 

27.8 

264. 

13160 

999.9 

999.9 

999.9 

27.3 

262. 

13430 

999.9 

999.9 

999.9 

25. 2 

263. 

13710 

150.0 

-56.3 

.000 

999.9 

999.9 

139  30 

999.9 

999.9 

999.9 

24.  7 

265. 

13990 

999.9 

999.9 

999.9 

24.2 

264. 

14280 

999.9 

999.9 

999.9 

25.7 

263. 

14570 

130.0 

-59.6 

.000 

999.9 

999.9 

148  30 

999.9 

999.9 

999.9 

24.7 

262. 

14860 

999.9 

999.9 

999.9 

23.2 

264. 

15160 

999.9 

999.9 

999.9 

22.7 

264. 

15470 

1 17.0 

-58.  0 

.000 

999.9 

999.9 

15500 

999.9 

999.9 

999.9 

20 . 6 

264. 

15730 

108.0 

-59.4 

.000 

19.0 

265. 

16000 

999.9 

999.9 

999.9 

17.  0 

259. 

16300 

100.0 

-58.3 

.000 

999.9 

999.9 

16480 

999.9 

999.9 

999.9 

16.  0 

266. 

16590 

999.9 

999.9 

999.9 

14.4 

264. 

16870 

999.9 

999.9 

999.9 

13.9 

264. 

17150 

999.9 

599.9 

999.9 

13.4 

269. 

17430 

999.9 

999.9 

999.9 

12.  4 

270. 

17710 

999.9 

999.9 

999.9 

10.3 

266. 

17980 

999.9 

999.9 

999.9 

9.8 

266. 

18260 

999.9 

999.9 

999.9 

9.3 

266. 

18540 

70.0 

-61.0 

.000 

999. 9 

999.9 

18710 

999.9 

999.9 

999.9 

8.2 

257. 

18820 

999.9 

999.9 

999.9 

8.2 

261. 

19100 

995.9 

999.9 

999.9 

5.7 

266. 

19380 

62.0 

-58.2 

.000 

999.9 

999.9 

19470 

60.0 

-58.5 

.000 

4.6 

262. 

19670 

999.9 

999.9 

999.9 

4.  1 

266. 

19960 

995.9 

99  9.9 

999.9 

4.6 

290. 

20240 

999.9 

999.9 

999.9 

6.7 

301. 

20520 

999.9 

999.9 

999.9 

4.6 

317. 

20800 
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V 

V 


50.0 

-53.9 

.000 

999.9 

999.9 

999.9 

999.9 

999.9 

999.9 

999.9 

999.9 

999.9 

42.0 

-54.3 

.000 

999.9 

999.  9 

999.9 

39.0 

-52.0 

.000 

999.9 

999.  9 

999.9 

36.0 

-52.6 

.000 

999.9 

999.9 

999.9 

999.9 

999.9 

999.9 

999.9 

999.9 

999.9 

30.0 

-48.  2 

.000 

999.9 

999.9 

999.9 

28.0 

-47.0 

.000 

999.9 

999.9 

20830 

4.1 

333. 

21120 

2.  1 

356. 

21440 

1.0 

034. 

21760 

999.9 

999.9 

21950 

1.0 

064. 

22080 

1.0 

267. 

22430 

0.5 

347. 

22770 

999.  9 

999.9 

22940 

0.5 

007. 

231 10 

0.  5 

238. 

23450 

0.5 

229. 

23790, 

1.5 

229. 

24130 

1.5 

207. 

24430, 

1.  5 

195. 

24590 
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APPENDIX  B 

SUMMARY  OF  SOBEOOTINES 

Subroutine  name;  IN  TIGS 

Description;  INTLGS  integrates  the  functions  ZL,  FL ,  CNSQI, 
CDSQI,  CQSQL  and  EPS  INI  over  the  slab 

Inputs;  P,  T,  Q,  STEITY,  SHEAR,  DQDZ,  RI,  THKSIB,  ZIMIN, 
ZIMAX,  Z1MEAN,  SGMSH1,  SGMST1 

Output;  PI,  PF,  PCNSQ,  PCDSQ,  PCQSQ,  PEPSIN 

Calls;  SGMSHI,  FNCTN1,  GSSIGN 

Called  by;  main  program 


Subroutine  name;  DATPRO 

Description;  DATPRO  interpolates  for  missing  rawinsonde 
data  and  calculates  the  weighted  average  and  gradient  values 

Inputs;  PRESS,  TEMP,  HUMID,  KS,  WD ,  Z,  N,  J,  RON 

Output;  PARY,  TARY,  QARY,  STBARY,  SHEARY,  DQDZRY,  RIARY, 
THKARY,  N 

Calls;  INTER 

Called  by;  main  program 


Subroutine  name;  INTER 

Description;  INTER  performs  a  semi-log  interpolation  of  the 
input  variables 

Inputs;  A,  B,  C,  D,  NUMBER 


45 


Output;  A NS 


Calls;  none 
Called  by;  DATPRO 


Subroutine  name;  FNC1NL 

Description;  FNCTN1  computes  the  value  of  the  turbulence 
parameters  (F,  Cn2,  and  EPSILON)  for  a  slab 

Inputs;  P,  T,  Q,  STBLTY,  SHEAR ,  DQDZ,  El,  THKSLB,  ZL, 
SGHSHR,  SGMSTB ,  ZL  ME  AN 

Output;  FL,  CNSQL,  CDSQL,  CQSQL,  EPSLNL 
Calls;  BVRT 
Called  by;  INTLGS 


Subroutine  name;  BVBT 

Description;  BVRT  controls  the  computation  of  the  turbu¬ 
lence  parameters  using  a  bivariate  model 

Inputs;  P,  T,  Q,  STBLTY,  SHEAR,  DQDZ,  RI,  ZL,  SGMSHR , 
SGMSTB 

Output;  F,  CNSQ,  CDSQ,  CQSQ,  EPSILN 
Calls;  INTGSH 
Called  by;  FNCI NL 


Subroutine  name;  INTGSH 


Description;  INTGSH  integrates  the  input  functions  over 
SHEAS  using  Gaussian  quadrature 


Inputs; 

SHRSLB, 

STBSLB, 

ZL,  SGHSHR,  SGMSTB 

Output; 

F,  CN2 , 

EPS 

Calls ; 

GSSLGN, 

FONSHE, 

INTGS2,  GSSLGR 

Called 

by;  BVRT 

Subroutine  name;  FUNSHR 

Description;  F  ON  SHE  computes  a  function,  FTAO,  from  a 
Rice-Nakagami  distribution  of  SHEAR'S  and  a  normal  distribu¬ 
tion  of  STABILITY 

Inputs;  TAU,  TAUSIB 

Output;  FENTAO 

Calls;  HBSS 10 

Called  by;  INTGSH 


Subroutine  name;  INTGS2 

Description;  INTGS2  integrates  the  normal  distribution 
function  together  with  (STABILITY)2  times  the  normal  distri¬ 
bution  function  over  STABILITY  using  Gaussian  quadrature 

Inputs;  STERIC,  STBSLB,  SGMSTB 

Output;  CDF3T3,  CN2STB 

Calls;  STBINT,  GSSLGR 


Called  by; 


INTGSH 


Subroutine  name;  STBINT 

Description;  STBINT  integrates  two  input  functions  over 
STABILITY  using  Gaussian  quadrature 

Inputs;  ST3MAX,  STBHIN,  STBSLB,  SGHSTB,  N STB,  NUMBEB 
Output;  CDFSTB,  CN2STB 
Calls;  GSSLGN,  EXPNX 
Called  by;  INT3S2 

Subroutine  name;  GSSLGN 

Description;  GSSLGN  calculates  the  abscissas  and  weights 
for  Gaussian  quadrature  for  limits  from  -1  to  +1 

Inputs;  N 

Output;  X,W 

calls;  none 

Called  by;  IN TGSH ,  INTGLS 

Subroutine  name;  GSSLGfi 

Description;  GSSLGB  calculates  the  abscissas  and  weights 
for  Gaussian  quadrature  for  limits  from  zero  to  infinity 

Inputs;  N 

Output;  X,$J 

Calls;  none 

Called  by;  INTGSH,  INTG52 


Subroutine;  MBS  SI  0 

Description;  BBSS  10  computes  the  value  of  the  modified 
Bessel  function  of  the  first  kind  of  order  zero  using  the 
IMSL  routine  MMBSIO 

Inputs;  X 

Output;  ANSWER 

Calls;  MBSSIO  (IMSL  routine) 

Called  by;  FONSHB 


Function  name;  EXPND 

Description;  EXPND  computes  the  value  of  a  normalized  expo¬ 
nential  distribution  for  ZL 

Inputs;  ZL#  ZLMEAN 

Output;  EXPND 

Calls;  EXPNX 

Called  by;  not  in  use  presently,  before  modification,  it 
was  called  in  FNCTNL 


Function  name;  SHRSGM 

Description;  SHRSGM  computes  the  standard  deviation  of  the 
SHEAR  distribution 

Inputs;  SGMCCN ,  P,  T,  ALPHA,  ST3LTY 

Output;  SHRSGM 

Calls;  none 

Called  by;  main  program 


Function  name;  STBSGM 

Description;  STBSGM  computes  the  standard  deviation  of  the 
fluctuating  part  of  STABILITY 

Inputs;  SHBSGM,  STBLTY 

Output;  STBSGM 

Calls;  none 

Called  by;  main  program 


Function  name;  SGMSH1 

Description;  SGMSHL  computes  the  value  of  SGMSHL  given 
SGMSH1  and  ZL 

Inputs;  SGMSH 1 ,  ZL 

Output;  SGMSHL 

Calls;  none 

Called  by;  INTLGS 


Function  name;  PROB19 

Description;  PR0B19  is  intended  to  compute  the  probability 
function  using  polynomial  approximation 

Inputs;  X 

Output;  PE0B19 

Calls;  none 


50 


Called  by;  not  presently  in  use,  before  modification,  it 
was  called  in  INTGSH 


Function  name;  EXPNX 

Description;  EXPNX  checks  the  input  arguments  to  ensure 
they  are  within  the  operating  limits  of  the  computer.  If 
they  are  not,  it  revalues  them  so  that  they  are  within 
limits  to  avoid  underflows  and  overflows 

Inputs;  X 

Output;  EXPNX 

Calls;  none 

Called  by;  ST3INT,  EXPFSH,  EXPND 


APPENDIX  C 

5ASPLE  T0BB2  DATAOOT  FILE  (iROHS  =  1) 


SLAB#  HIDA1T  { M) 


1 

1645.00 

2 

1805.00 

3 

2065.00 

4 

2220.00 

5 

2380.00 

6 

2595.00 

7 

2885.00 

8 

3120.00 

9 

3350.00 

10 

3640.00 

1 1 

3825.00 

12 

4035.00 

13 

4185.00 

14 

4310.00 

15 

4425.00 

16 

4540.00 

17 

4650.00 

18 

4790.00 

19 

4930.00 

20 

5050.00 

21 

5230.00 

22 

536  5.00 

23 

5480.00 

24 

5620.00 

25 

5720.00 

26 

5865.00 

27 

5995.00 

28 

6135.00 

CN2 

0.  6 42 9 89 5 07 6D-  1  6 
0. 8525603886 D- 16 
0.  1  3087 13280D- 1  6 
0.  1840405935D-16 
0. 1675688171D- 1 6 
0.  9289 98 0295 D- 17 
0.  9  555537 126D- 1 7 
0. 1028519061D-16 
0.  1845121 879D- 1 7 
0-8  74  6728887  D- 17 
0.  2248982  453D- 1 6 
0. 7019746893D-1 6 
0.  167 6 746 73 6D- 1 5 
0.207  62 5 3 79 2D- 1 5 
0. 1 877315078D- 1 5 
0. 2996434361 D- 16 
0-  2 88 49 30 02 8D-  1  6 
0. 1875013105 D- 16 
0.  2  125610311D-16 
0.2529 149526D- 16 
0. 2375408888D- 1 6 
0.1080201618D-16 
0. 1088445482D- 1 6 
0. 8812439910D-1 7 
0.  9 48 6 5 86 7 8 9D-  1  7 
0.  1  078625 1 92D-1 6 
0.  1777040099D- 1 6 
0. 2783679743D- 1 6 
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29 

6320.00 

30 

6435.00 

31 

6500.00 

32 

6575.00 

33 

6690.00 

34 

6750.00 

35 

6870.00 

36 

7105.00 

37 

7300.00 

38 

7390.00 

39 

7460.00 

40 

7560.00 

41 

7705.00 

42 

7960.00 

43 

8275.00 

44 

8495.00 

45 

8640.00 

46 

8855.00 

47 

9030.00 

48 

9150.00 

49 

9325.00 

50 

9445.00 

51 

9490.00 

52 

9595.00 

53 

9745.00 

54 

9870.00 

55 

10075.00 

56 

10360.00 

57 

10575.00 

58 

10710.00 

59 

10885.00 

60 

11010.00 

61 

11150.00 

62 

11410.00 

63 

11665.00 

0.31561 49385D- 1 6 
0. 23737861 14D-1 6 
0.  224993702  ID-  1  6 
0.2985857  815D-1 8 
0.  1865800713D-1 6 
0. 1937612712D-16 
0. 1 052763 1 86D- 1 6 
0.1 490718353D-1 6 
0.  1813311 7 19D- 1 6 
0.2 892384  730D-1 6 
0.  27450  28  922D- 1 6 
0. 1685917667D-16 
0.  1 6461 31 759D- 1 6 
0.  1 687852523D-1 6 
0.  8895951 535D- 1 7 
0. 1 192759561D-1 6 
0.  1  121599474D-1 6 
0.9671650870D-17 
0. 1045068773D- 16 
0. 1359829831D-16 
0.3 27 9598721 D-1 6 
0. 3337614 124D- 16 
0.  1042557235D-  1  5 
0.6504480 124 D-1 6 
0. 4285484888D- 1 6 
0.4098830537 D-1 6 
0. 3 1 10753455D- 1 5 
0. 1411441 978D- 1 6 
0.  1  27  24  54  52 2D-  1  6 
0. 1333465354D-1 6 
0. 583 1460262D- 1 6 
0.5737078697 D-1 6 
0. 3580049158D- 16 
0.  1 080  384642D- 1 6 
0. 9533808 124D- 1 7 
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11920.00 

0.  1521933891D-16 

12075.00 

0. 1874242331D- 1 6 

12210.00 

0.1619 82 9 77  ID- 16 

12460.00 

0.2374495506D-16 

12740.00 

0.2598101284D-16 

13020.00 

0. 1549131 182D-16 

13295.00 

0. 1054737223D-16 

13570.00 

0.  15401 70  661D- 1 6 

13820.00 

0. 1251849086D-16 

13960.00 

0. 1213451 593D- 1 6 

14135.00 

0. 1147173704D-16 

14425.00 

0. 5974197139D-1 8 

14700.00 

0. 1379922626D-1 7 

14845.00 

0. 1774158 125D- 1 7 

15010.00 

0. 28833704 86D-17 

15315.00 

0.  1067950  577D- 1 6 

15485.00 

0.1 216896389D-1 6 

15615.00 

0.9765411 42 2D- 1 7 

15865.00 

0.  8095732044D-1  7 

16150.00 

0.  1 096  7630  98D- 1 6 

16390.00 

0. 9651295 96 8D- 17 

16535.00 

0.  9  20  7633  50 3D-  1  7 

16730.00 

0.8188833993D-17 

17010.00 

0. 6958337906D- 1 7 

17290.00 

0.7  09  0243  29  2D- 17 

17570.00 

0. 6583973759D- 1 7 

17845.00 

0.7582720669D-17 

18120.00 

0. 58361 68579D- 1 7 

18400.00 

0. 4927868104 D- 17 

18625.00 

0.  5 6808 6 2 26 8D-  1  7 

18765.00 

0. 5600901 504 D- 17 

18960.00 

0.  4850331 731D-1 7 

19240.00 

0. 7716729501 D- 17 

19425.00 

0.6594158104D-17 

19570.00 

0.4 02 3897304 D- 17 

99 

19815.00 

0.  373464229 9D- 1 7 

00 

20100.00 

0. 4 26 132 0896 D-1 7 

01 

20380.00 

0. 4530608595D- 1 7 

02 

20660.00 

0. 5746 32 4061 D-1 7 

03 

20815.00 

0.4799269  686D- 1 7 

04 

20975.00 

0. 4  589887089D-1 7 

05 

21280.00 

0.  4683289  602D- 1 7 

06 

21600.00 

0.2859715419 D-1 7 

07 

21855.00 

0.  255  1690 196D- 1 7 

08 

22015.00 

0. 30959744 97D- 17 

09 

2225  5.00 

0. 3  2294 1 1 536D- 1 7 

10 

22600.00 

0.2199815305D-17 

11 

22855.00 

0.2  34  2476  57  ID- 17 

12 

23025.00 

0. 2279193328 D- 17 

13 

23280.00 

0.  201 5 465 64 4D- 1 7 

14 

2362  0.00 

0. 2460930328D-1 7 

15 

23960.00 

0. 2368559736D- 1 7 

16 

24280.00 

0. 1883079643D-17 

17 

24510.00 

0.  1813421 2G4D- 1 7 

cwxmm 


APPEND!!  D 

CCHPAHISOH  OP  BESDLTS 


Table  4  illustrates  the  corresponding  values  of  Cn2  (z) 
obtained  from  the  NOAA  program  and  the  IBM  model  at  the 
specified  altitudes.  The  fourth  column  lists  the  percent 
difference  in  the  IBM  value  relative  to  the  NOAA  value.  The 
#EONS  was  set  at  1  which  resulted  in  the  MIDALT  values  being 
the  same  for  both. 


TABLE  4 

Comparison  of  Program  Output  with  *BUNS  =  1 


MIDAIT  (m) 

NOAA  Cn2 

MODEL  Cn2 

%  DIFFERENCE 

1810.00 

0.7908E-16 

0. 8526E-16 

7.8  IS 

2070.00 

0. 1029E-16 

0. 1309E-16 

27.18% 

2220.00 

0.518  9E-  2  0 

0. 1840E-16 

354574.06% 

2380.00 

0.5120E-20 

0.  1676E- 1 6 

327182.81% 

2600.00 

0.108  IE- 16 

0.  9290E-1 7 

14.06% 

2880.00 

0.9  127E- 17 

0.  9556E-17 

4.70% 

3120.00 

0.140  9E- 1  8 

0. 1029E-16 

7199.64* 

3350.00 

0. 2106 E- 18 

0.  1845E-17 

776.13% 

3640.00 

0.379  3E- 1  8 

0. 8747E-1 7 

2206.02% 

3830.00 

0. 2510  E- 16 

0. 2249E-16 

o 

9 

o 

4040.00 

0.566  0E- 16 

0. 7020E-16 

24.02% 

4190.00 

0. 1529E-15 

0. 1677E-15 

9.66% 

4310.00 

0.208  0E-  1  5 

0. 2076E-1 5 

0.18% 

4420.00 

0.207  IE-15 

0. 1877E-15 

9.35% 

4540. 00 

0.2880E-16 

0. 2996E- 1 6 

4.04% 

4650.00 

0.3446E-16 

0.  2885E-16 

16.28% 

4790.00 

0. 1726E-16 

0. 1875E-16 

8.63% 

4930.00 

0.194  2E- 16 

0.2126E-16 

9.45% 

5050.00 

0.2452E-16 

0.  2529E-1 6 

3.15% 

5230.00 

0.22 8 8E- 16 

0.2375E-16 

3.82% 

5360.00 

0. 1090E-16 

0. 1080E-16 

0.90% 

5480.00 

0. 1029E-16 

0. 1 088E-1 6 

5.78% 

5620.00 

0.8983E-17 

0. 8812E-1  7 

1.90% 

5720.00 

0.917  7E- 1 7 

0. 9487E-1 7 

3.37% 

5870.00 

0. 106  3E-16 

0. 1079E-16 

1  .47% 

6000.00 

0.16 19E-16 

0. 1777E-16 

9.76% 

6140.00 

0.270  5E- 16 

0. 2784E-16 

2.91% 

6320.00 

0.2939E-16 

0.  3156E-16 

7.39% 

6430.00 

0.3485E-20 

0. 2374E-1 6 

681043.31% 

6500.00 

0.346  7E-20 

0.  2250E-16 

648857.37% 

6570. 00 

0.1168E-20 

0.  2986E-1 8 

25463.34% 

6690.00 

0.  1812E-16 

0.  1866E-1 6 

2.97 % 

6750.00 

0 . 188  7E-  16 

0.  1  938E-16 

2.68% 

6870.00 

0. 1089  E- 16 

0.  1053S-16 

3.3  3% 

7100.00 

0.146  2E- 1 6 

0. 1491E-16 

1  .96% 

7300.00 

0. 1720E-16 

0.  1813E-16 

5.43% 

7390.00 

0-197  4E- 16 

0. 2  892E-1 6 

46.52% 

7460.00 

0.18 12E-16 

0.  2745E-16 

51.49% 

7560.00 

0 . 1  6 1 5E-  1 6 

0. 1686E-16 

4.39% 

7700.00 

0.8993E-17 

0. 1646E-16 

83.05% 

7960.00 

0.923  IE-  17 

0.  1688E-16 

82.85% 

8280.00 

0. 6828  E- 17 

0.  8896E-17 

30.29% 

8500.00 

0.1 055E-16 

0. 1 193E-16 

13.06% 

8640.00 

0. 8701 E- 17 

0. 1122E-16 

28.90% 

3850.00 

0.5823E-17 

0. 9  672E- 1  7 

66.09% 

9030.00 

0.592SE-17 

0.  1045E-16 

76.26% 

9150.00 

0.133  5E-  16 

0. 1360E-16 

1.86% 

9330.00 

0 . 355  9E-16 

0. 3280E-16 

7.85% 

9440.00 

0.3  479E- 16 

0. 3338E-1 6 

4.06% 

9490.00 

0.9  17  5E-  16 

0.  1043E-1  5 

13.63% 

9600.00 

0.581 2E-  1 6 

0. 6504E-1 6 

11.91% 
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17010.00 

0.5884E-17 

0. 6958E-17 

18.26% 

17290.00  . 

0.6150E-17 

0. 7090E- 17 

15.  29% 

17570.00 

0.5682E- 17 

0. 6584E-17 

15.87% 

17840.00 

0.6798E-17 

0. 7583E- 17 

11.54% 

18120.00 

0.4908E-17 

0. 5836E- 17 

18.91% 

18400.00 

0.466  IE-17 

0.4928E-17 

5.73% 

18620.00 

0.5210E-17 

0. 568  IE-  1 7 

9.04% 

18760.00 

0.5  867E- 17 

0. 560 1 E- 17 

4.54% 

18960.00 

0.5052E-17 

0. 4850 E-  17 

3.  99% 

19240.00 

0.640  IE-  17 

0. 7717E-17 

20.  56% 

19420.00 

0.4795E-17 

0. 6594E-17 

37.52% 

19570.00 

0.383  8E- 1 7 

0. 4024E-17 

4.  84% 

19810.00 

0.436  2E-17 

0 . 3735E-  17 

14.38% 

20100.00 

0.473  5E-  1 7 

0. 4261E-17 

o 

o 

t 

o 

20380.00 

0. 4838E-17 

0.  45 3 IE- 17 

6.  35% 

20660.00 

0.4805E-17 

0. 5746E-  17 

19.59% 

20820.00 

0. 3779 E- 17 

0.4799E-17 

27. 00% 

20970.00 

0.312  2E- 1 7 

0.4590E-17 

47. 02% 

21280.00 

0.3422E-17 

0.  4683E-  17 

36.  86% 

21600.00 

0.282  4E-17 

0.2860E-17 

1.26% 

21850.00 

0 . 2532E-17 

0.  2552E- 1 7 

0.78% 

22020.00 

0.302  IE-17 

0. 3096E-17 

2.  48% 

22250.00 

0.3172E-17 

0. 3229E-17 

1.81% 

22600.00 

0.2136E-17 

0. 2200E-  17 

2.  99% 

22860.00 

0. 1972E-17 

0.  2342E-  17 

18.  79% 

23030.00 

0.241 5E-  17 

0. 2279E-17 

5.  62% 

23280.00 

0.2358E-17 

0. 2015E-17 

14.  53% 

23620.00 

0.216  9E- 1 7 

0 . 246 1 E-  17 

13.46% 

23960.00 

0. 2084  E- 17 

0.  2369E-17 

13.65% 

24280.00 

0. 1840E-17 

0. 1883E-17 

2.34% 

Table  5  illustrates  the  corresponding  values  of  Cn2 (z) 
obtained  from  the  NOAA  program  and  the  IBM  model  at  the 
specified  altitudes.  The  fifth  column  lists  the  percent 
difference  between  the  IBM  value  and  the  interpolated  (semi¬ 
log)  value  from  the  NOAA  data.  Note  that  the  different 
interpolation  method  (i.e.  IRONS  =  2)  produced  different 
values  for  the  MIDALT*s  thus  making  direct  comparison 
impossible. 
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